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ABSTRACT 


A  data  acquisition  and  reduction  system  has  been  created 
for  aircraft  materials  fatigue  testing.  The  system  uses  an 
HP-9835  Desktop  Calculator,  an  HP-2240A  Measurement  and 
Control  Processor  and  a  Materials  Testing  System  loading 
machine.  Three  different  computer  programs  on  the  HP-9835 
are  used  to  analyze  material  properties,  simulate  inflight 
fatigue  loading  and  compute  fatigue  damage  at  a  stress  con¬ 
centration.  The  flight  loads  are  selected  from  Mil  Spec  8866 
and  applied  in  a  random  order.  The  fatigue  damage  at  a  stress 
concentration  is  calculated  from  the  applied  local  stresses 
using  Miner's  Law. 
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I .  INTRODUCTION 


Today's  high  cost  of  aircraft  dictates  that  the  lifetimes 
of  existing  and  new  aircraft  be  extended  as  long  as  possible. 
The  data  acquisition  system  and  computer  programs  described  in 
this  thesis  will  permit  more  realistic  fatigue  testing  of 
aircraft  materials;  thereby  enabling  the  engineer  to  predict 
and  design  more  accurately  for  fatigue  life. 

The  programs  are  set  up  to  start  with  a  standard  tension 
test  to  determine  the  modulus  of  elasticity  and  the  Ramberg- 
Osgood  coefficients  to  model  the  plastic  region.  The  second 
program  is  for  flight  load  simulation  of  a  specific  number 
of  flight  hours  to  be  applied  to  a  specimen  with  a  given 
stress  concentration.  This  program  was  designed  to  apply 
loads  that  would  be  experienced  by  fighter  aircraft,  the 
loads  being  randomly  ordered  after  being  selected  from  the 
Military  Specification  8866  1000  flight  hour  spectrum.  The 
load  data  from  this  test  is  stored  in  sequence  on  magnetic 
tape,  and  local  strain  information  is  gathered  at  the  point 
of  stress  concentration  via  pictures  of  photo-elastic 
patterns  taken  during  the  test.  The  last  program  assesses 
the  damage  incurred  by  the  specimen  during  the  flight 
simulation.  This  program  is  specialized  to  predict  damage 
for  specimens  fabricated  from  7075-T6  Aluminum. 
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II.  CONCEPTS 


A.  STRESS  AND  STRAIN  CALCULATIONS 

Lindsey  has  developed  [Ref.  1]  a  series  of  equations 
for  calculating  local  stresses  at  stress  concentration 
sites,  which  use  a  nonlinear  stress-strain  law  and  account 
for  residual  stresses  and  strains.  The  analysis  assumes 
that  the  notch  tip  behaves  in  many  respects  like  a  uni¬ 
axial  tension  specimen  with  a  Ramberg-Osgood  description. 


e 


3(|)n 


(1) 


Loading  the  specimen  of  Figure  1  to  a  nominal  stress  of 
S  psi,  the  local  stress  can  be  found  by  solving 


KtS  =  o[l  +  3(|)n'1]  (2) 

where 

Kfc  =  Elastic  Stress  Concentration  Factor 
During  the  unloading  portion  of  the  cycle,  Lindsey  showed, 
via  energy  considerations,  that  the  residual  stress  was 
given  by 

to2  ♦  «|,nV 

where  or  is  the  local  stress  at  the  notch  at  the  peak  of  the 
loading  cycle  as  found  from  Equation  (2)  above.  This  stress 
will  be  compressive  and  will  be  non-zero  only  if  a  exceeds 
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ujure 


the  proportional  limit  by  a  significant  amount.  The 
residual  strain  is  positive  and  given  by 


= 


_R 

E 


3(|)n 


(4) 


These  values  of  stress  and  strain  become  the  starting  points 
for  the  next  cycle,  and  the  R  fatigue  parameter  is  calculated 
taking  the  residual  stress  into  account.  (The  R  value  is 
the  ratio  of  the  minimum  stress  to  the  maximum  stress.)  The 
elastic  stress  in  subsequent  cycles  is  given  by 


c  =  +  Krs  -  EeR  [1  ♦  ,  1 


(5) 


Yield 

At  the  next  yield,  a  new  residual  stress  and  strain  must  be 
calculated  from  Equations  (3)  and  (4) . 


B.  MILITARY  SPECIFICATIONS  FOR  FATIGUE 

In  order  to  analyze  materials  and  structures  for  fatigue, 
it  is  necessary  to  have  some  knowledge  of  the  types  of  loads 
and  the  number  of  loading  cycles  that  the  structure  will  be 
subjected  to.  For  this  purpose  Military  Specification  8866 
Flight  Loading  Spectrum  A  has  been  chosen  (Table  I) .  This 
spectrum  was  chosen  by  the  military  for  loadings  that  were 
expected  on  fighter-type  aircraft  for  each  1000  hours  of 
flight  time.  The  military  expects  these  aircraft  to  have 
a  lifetime  of  4000-6000  hours  under  this  type  of  loading. 


TABLE  I 


FLIGHT  LOAD  SPECTRUM  A 


%  Limit  Load 

Number  of  Cycles 

35 

17,000 

45 

9,500 

55 

6,500 

65 

4,500 

75 

2,500 

85 

1,360 

95 

440 

105 

150 

115 

40 

125 

16 

C.  STANDARD  USED  FOR  MATERIAL  PROPERTIES 

The  reference  material  specifications  for  fatigue  life 
were  taken  from  NACA  TN  2324  [Ref.  2].  Testing  was  done  on 
sheet  75S-T6  Aluminum  (A1  7075-T6)  for  various  R  values. 

The  results  of  these  tests  were  plotted  in  Figure  2. 
Regression  equations  were  formulated  using  the  BMDP 
Biomedical  Computer  Programs  [Ref.  3]  by  inputing  the  data 
points  in  Table  II  that  were  selected  off  Figure  2. 

The  best  regression  formula  (Table  III)  was  found  by 
using  the  actual  R  value  and  the  stress  divided  by  10,000. 
The  output  value  was  the  logarithm  of  the  number  of  cycles 
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S  (KSI) 


*  1 


TABLE  III 


REGRESSION 

COEFFICIENTS 

Term 

Coefficient 

INTERCEPT 

12.6452 

S 

-1.92662 

s4 

2.81098E-3 

s8 

-3. 10691E-7 

R2 

-12.8099 

R6 

212. 476 

R8 

-86.7650 

SR 

3.68800 

s3r 

-0.112720 

s5r 

1.04762E-3 

sV 

3 . 39637E-6 

s3r4 

-3 . 50885E-2 

s4r6 

-1. 61827E-2 

SR9 

-34.4642 

S  =  Stress/10000 

R  =  Minimum  Stress/Maximum  Stress 

to  failure  (N) .  Thus,  the  logarithm  of  N  is  the  sum  of  the 
terms  shown  in  Table  III.  This  gave  a  multiple  correlation 
coefficient  of  0.96167.  Table  IV  shows  the  correlation  of 
the  regression  output  to  the  input  values  of  Table  II. 

16 
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Although  the  regression  predicted  output  does  not  look  like 
a  close  correlation  to  the  observed  data,  the  error  of  the 
regression  equation  falls  within  the  scatter  band  for  the 
original  data  [Ref.  2]. 
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III.  DESCRIPTION  OF  ACQUISITION  SYSTEM 


The  data  acquisition  system  consists  of  an  HP-9835 
Desktop  Computer  (Controller) ,  an  HP-2240A  Measurement  and 
Control  Processor  (Processor) ,  and  a  Materials  Testing 
Machine  (MTS) .  The  Controller  (Fig.  3)  is  used  to  generate 
the  random  loads,  send  the  loads  to  the  MTS  through  the 
Processor  and  manipulate  the  strain  data  that  it  receives 
back  through  the  Processor.  The  Processor  (Fig.  4)  receives 
digital  commands  from  the  Controller  to  output  specified 
analog  values  or  to  read  analog  values  and  send  them  in 
digital  form  to  the  Controller.  The  MTS  (Fig.  5)  is  a 
hydraulic  testing  machine  that  applies  loads  specified  by 
voltage  signals  to  a  test  specimen. 

The  resident  language  in  the  Controller  is  BASIC  in 
which  all  of  the  programs  are  written.  Programs  and  data  are 
stored  on  a  built-in  magnetic  tape  mass  storage  device.  The 
Controller  is  linked  to  the  Processor  (or  any  other  peripherals) 
through  the  HP-IB  Interface  Bus.  This  Bus  has  its  own  Read 
Only  Memory  (ROM)  and  address  associated  with  the  ROM.  Each 
peripheral  on  the  Bus  also  has  its  own  HP-IB  address.  The 
Processor  is  addressed  by  the  Controller  through  the  HP-IB 
by  use  of  the  ENTER  and  OUTPUT  commands,  which  include  the 
device  address  and  the  data  to  be  sent  or  requested. 
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Figure  5.  Materials  Testing  System  (MTS) 


The  Processor  has  four  function  card  slots.  For  this 
data  acquisition  system  one  D/A  card  and  one  A/D  card  were 
all  that  were  necessary.  The  D/A  card  can  handle  up  to  four 
output  channels,  one  is  used  to  drive  the  MTS  and  another  is 
used  to  drive  a  relay,  which  operates  an  automatic  camera 
shutter.  The  A/D  card  can  handle  up  to  sixteen  bi-polar 
channels.  In  this  instance  eight  were  used  by  the  acquisition 
program  for  reading  strain  gages,  and/or  load  cells.  Commands 
are  sent  to  the  Processor  in  standard  ASCII  format,  and  an 
error  code  is  sent  back  to  the  Controller. 

The  MTS  is  a  100,000  lb.,  closed  loop,  feedback  controlled 
Materials  Test  System  with  a  hydraulic  actuator.  The  MTS  can 
be  used  in  either  load,  strain  or  displacement  control  modes. 
It  can  be  used  by  itself  with  its  own  controller  or  its  own 
function  generator.  In  this  case  the  output  from  the  D/A 
card  was  input  to  the  MTS  in  a  load  control  mode.  The  signal 
from  the  Processor  was  output  to  the  hydraulic  system  and  the 
MTS  checks  itself  with  a  load  cell  signal  until  the  load  is 
matched. 
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IV.  DESCRIPTION  OF  PROGRAMS 


There  are  three  major  programs  written  for  the  controller 
(in  BASIC)  for  material  properties  and  fatigue  testing.  The 
first  is  entitled  AXIAL  and  is  used  to  record  data  from  a 
strain  gaged  specimen  in  an  axial  tension  test  to  determine 
the  Modulus  of  Elasticity/  Poisson's  Ratio  and  the  Ramberg- 
Osgood  coefficients  for  the  plastic  region.  FLTTST,  the 
second  program,  generates  random  ordered  loads  and  outputs 
these  loads  through  the  PROCESSOR  to  the  MTS  System,  which 
applies  the  loads  to  the  specimen.  The  loads  and  the  strains 
generated  by  this  program  are  stored  on  tape  for  use  later  in 
data  reduction.  The  last  program,  DAMAGE,  takes  the  data 
generated  by  FLTTST  and  calculates  the  fatigue  damage  of  the 
specimen. 

A.  PROGRAM  AXIAL 

AXIAL  does  not  control  the  MTS  System  and  so  can  be  used 
with  any  system  possessing  voltage  outputs  for  load  and 
strain  to  record  data.  AXIAL  records  data  while  the  testing 
machine  performs  an  axial  tension  pull  on  a  test  specimen. 

The  program  is  set  up  to  read  the  load  and  both  the  longi¬ 
tudinal  and  transverse  strains  at  equal  intervals  of  time 
while  the  testing  machine  is  used  in  a  constant  strain  rate 
mode.  The  program  can  also  be  run  using  only  the  longitudi¬ 
nal  strain  gage  in  which  case  Poisson's  ratio  is  not  computed. 
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The  loads  ax~  converted  to  stresses  so  that  the  program  can 
work  directly  with  the  stress-strain  relationships.  After 
all  the  data  has  been  taken,  the  program  asks  the  user  to 
specify  the  load  at  which  the  proportional  limit  seems  to 
have  been  reached.  This  cutoff  is  entered  so  that  a  good 
linear  fit  can  be  made  to  find  the  Modulus  of  Elasticity  (E) . 
After  the  linear  fit  to  obtain  E,  the  value  is  output  along 
with  the  slope  intercept  B  and  the  correlation  factor. 

The  program  then  asks  the  user  if  a  different  fit  is  to  be 
tried  using  a  new  proportional  limit.  If  not,  the  program 
continues  by  computing  Poisson's  ratio  and  finds  the  Ramberg- 
Osgood  coefficients  using  the  Modulus  of  Elasticity  previous¬ 
ly  found. 

The  Modulus  of  Elasticity  is  the  most  critical  property 
found  in  AXIAL  not  only  for  its  own  value  in  calculating 
stresses  and  strains  but  also  because  the  Ramberg-Osgood 
coefficients  are  based  on  E.  For  this  reason  a  good  linear 
fit  is  obtained  by  a  routine  that  evaluates  the  perpendicular 
distances  of  the  data  points  from  the  estimated  fit.  The 
E  that  has  be^n  calculated  is  output  to  the  user  along  with 
the  correlation  factor  so  that  a  judgment  can  be  made  as  to 
whether  the  estimated  proportional  limit  should  be  changed 
in  order  to  try  for  a  better  fit. 

The  equations  governing  the  linear  elastic  portion  of  the 
stress-strain  curve  are: 

7^  *  As^  +  B  (6) 

7  ^  £  2  ( 7 ) 
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where  A  is  the  modulus  of  elasticity  (E) .  B  and  D  are 
normally  neglectibly  small.  Poisson's  ratio  is  given  by 


The  plastic  portion  of  the  stress-strain  curve  can  be 
approximated  by  the  following  equation: 


e  -  |  +  3(f)n  (1) 

Beta  and  n  can  be  found  by  performing  a  logarithmic  curve 
fit  of  the  data  above  the  proportional  limit,  but  the 
Ramberg- Osgood  method  gives  as  good  a  fit  with  fewer  calcu¬ 
lations.  The  Ramberg-Osgood  relations  are: 


n 


i  +  .  ipa  - 

log  V7/a>85) 


(9) 


3  =  4(-^-)n_1  (10) 

7  J.7 

where , 

o  -j  -  stress  at  the  intersection  of  a  line  of  slope  . 7E 
and  the  stress-strain  curve 
a  QC = similarly  computed  with  a  slope  of  .85E. 

•  O  J 

Inputs  for  AXIAL  are: 

Maximum  Strain  for  the  test 
Strain  Rate  to  be  used 

Number  of  Data  Points  desired  (1000  max) 

Number  of  Strains  to  be  recorded  (1/2) 

Cross  Sectional  Area  of  specimen 
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Load  Scale  Factor  (pounds/millivolt) 

Strain  Scale  Factor  (microstrain/millivolt) 

After  the  data  run: 

Maximum  Load  of  Elastic  Range  (Proportional  Limit) 

B.  PROGRAM  FLTTST 

The  most  critical  program  is  FLTTST  since  it  both 
controls  the  MTS  and  takes  data.  The  key  to  FLTTST  is  the 
random  selection  of  the  load  to  be  applied  to  the  specimen. 
Using  the  Mil  Spec  8866  Flight  Loading  Spectrum  A  for 
fatigue,  percentages  were  calculated  for  the  number  of  times 
each  load  is  to  be  applied  (Table  V) .  These  percentages 
were  then  totaled  and  left  in  decimal  form  as  numbers 
between  0  and  1.  Then,  using  the  Controller's  random 
number  generator,  the  load  to  be  applied  is  selected  through 
a  series  of  IF  statements.  For  example,  to  select  a  load 
of  95%  of  limit  load,  a  random  number  above  0.98462125  and 
below  0.99509594  has  to  be  selected.  If  after  going  down 
through  all  of  the  IF  statements  a  load  is  not  selected,  it 
would  mean  that  the  random  number  must  be  above  0.99961910, 
then  125  limit  load  has  been  selected,  the  MTS  load  is 
incremented  in  finite  load  steps  that  are  user  selected. 

When  the  desired  load  is  reached,  the  Controller  than 
requests  that  the  Processor  read  the  data  channels.  FLTTST 
then  stores  these  data  points  on  tape  for  future  reduction. 
As  each  load  to  be  output  is  selected,  a  count  is  kept  of 
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that  load  so  that  the  user  can  have  an  immediate  output  of 
the  load  distribution  at  the  end  of  the  run. 

TABLE  V 

RANDOM  NUMBER  FOR  SELECTING  MANEUVER  LOAD 


Limit  Load 

#  Cycles 

Random  Number 

35 

17,000 

0.40470409 

45 

9,500 

0.63086226 

55 

6,500 

0.78560206 

65 

4,500 

0.89272961 

75 

2,500 

0.95224492 

85 

1,360 

0.98462125 

95 

440 

0.99509594 

105 

150 

0.99866686 

115 

40 

0.99961910 

Inputs  for  FLTTST  are: 

Number  of  Hours  to  be  simulated 
Number  of  Data  Readings  (maximum  of  8) 

Load  at  One  G  (lbs) 

Limit  Load  (lbs) 

Load  Scale  of  MTS  (10,000;20,000;50,000;100,000) 
Load  Rate  (lbs/sec) ,  and 

Data  Reading  Scale  Factor  (microstrain  or  pounds/ 
millivolt) 
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There  is  another  program  named  RDDATA  that  may  sometimes 
be  needed  after  running  FLTTST.  If  data  from  more  than  one 
strain  gage  or  load  cell  were  recorded  on  the  data  file 
from  FLTTST,  it  is  necessary  to  first  run  RDDATA  to  separate 
each  data  file  under  its  own  heading  before  proceeding  with 
program  DAMAGE.  This  was  done  to  conserve  time  between  data 
points  taken  and  to  extend  the  life  of  data  tapes.  It  also 
conforms  with  program  DAMAGE  since  it  can  only  compute 
damage  on  one  set  of  data  at  a  time. 

C.  PROGRAM  DAMAGE 

The  purpose  of  this  program  is  to  take  the  loading  that 
was  saved  on  tape  with  the  FLTTST  program  and  calculate  the 
fatigue  damage  accumulated  at  the  stress  concentration  site 
during  the  flight  hours  simulated.  The  program  reads  the 
recorded  far-field  loads  in  the  order  that  they  were  applied 
to  the  specimen,  calculates  local  stresses  at  the  stress 
riser  and  then  uses  Miner's  Law  for  calculating  the  cumula¬ 
tive  damage.  Once  the  load  creates  a  stress  in  the  critical 
area  greater  than  the  proportional  limit,  local  stresses  are 
calculated  using  a  nonlinear  stress-strain  law  and  calcula¬ 
tions  for  residual  stresses  and  strains  are  made. 

The  maneuver  load  and  the  one  G  load  stored  immediately 
after  it  are  read  in  together.  The  maneuver  load  is  then 
converted  to  nominal  or  far-field  stresses  and  local  stresses 
using  the  following  relations: 

3  *5  'll> 
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a  =  \S 


(12) 


where , 

S  =  nominal  stress 
L  =  applied  load 

A  =  minimum  cross  sectional  area  of  specimen 

0  =  local  stress 

Kt=  stress  concentration  factor 

The  local  maneuver  stress  is  checked  against  the  proportion¬ 
al  limit  or  the  last  largest  plastic  stress  to  see  if  it 
would  cause  plastic  deformation.  If  plastic  deformation  is 
produced,  then  an  iteration  process  must  be  performed  to 
find  the  local  stress,  using  the  relation 


KtS  -  a[l  +  6(|)n>'1]  (2) 

The  routine  to  do  this  uses  a  modified  false  position  method 
[Ref.  4].  With  the  local  stress  determined,  residual 
stresses  and  strains  are  then  calculated,  using 


p  4 

~  0  -  [0 


2ngE< 

n+1 


(I) 


n+1. 


(3) 


R 


3(|)n 


(4) 


If  no  new  plastic  deformations  are  caused  but  the  plastic 
region  has  been  entered  before,  then  local  stresses  are  cal¬ 
culated  using 


0 


7R  +  KtS 


Ecr[!  + 


SLast 

Yield 


(5) 
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where  ST  .  is  the  largest  nominal  stress  that  has  occurred 
Last 

Yield 

which  produced  local  stresses  into  the  plastic  region. 

Damage  is  then  calculated  using  the  current  local  stress 
and  the  ratio  of  the  current  stress  to  the  stress  at  one  G 
loading  (R  value) . 

Using  Miner's  Law: 


where,  N  *  number  of  cycles  to  failure  and  can  be  found  using 
regression  formulas  that  fit  stress-fatigue  data  for  7075-T6 
aluminum.  The  damage  caused  by  this  cycle  is  then  added  to 
the  damage  caused  by  the  other  cycles  computed  one  load  at 
a  time  until  all  the  data  has  been  read. 

Inputs  to  DAMAGE  are: 

Minimum  Cross  Sectional  Area  (sq.  in.) 

Modulus  of  Elasticity  (psi) 

Ramberg-Osgood  Coefficients,  and 
Proportional  Limit  (psi) 


31 


V.  SAMPLE  TEST 


Several  tests  have  been  developed  to  evaluate  both  the 
programs  and  the  program-acquisition  system  combined.  AXIAL 
was  checked  by  inputting  stress  and  strain  pairs  taken  from 
an  existing  Ramberg-Osgood  representation  of  7075-T6,  and 
the  results  of  the  two  were  compared.  FLTTST  was  checked 
by  taking  the  output  from  the  Processor  (normally  used  to 
drive  the  MTS) ,  routing  the  output  from  the  output  port  to 
the  input  port  and  reading  the  output  value  on  the  data 
input  port.  Thus,  the  input  port  was  getting  input  just  as 
it  normally  would  from  a  load  cell.  By  doing  this  the  actual 
voltages  could  be  read  directly  from  the  output  port  with 
a  voltage  meter,  and  the  input  results  could  be  treated  as 
actual  loads  to  check  Program  DAMAGE.  This  was  checked  by 
hand  for  both  load  spectrum  and  load  accuracy.  DAMAGE  was 
checked  by  hand  with  sample  runs  through  FLTTST  as  described 
above.  Two  full  system  tests  were  then  run  with  the  MTS, 
including  photo-elastic  pictures  of  the  stress  concentration 
area  and  strip  chart  recordings  of  the  input  and  output  loads. 

A.  PROGRAM  AXIAL 

The  first  test  of  AXIAL  was  performed  by  using  the 
Ramberg-Osgood  coefficients  calculated  by  Kaiser  [Ref.  5]  to 
generate  stress-strain  pairs  for  manual  input  into  the 


32 


program.  The  result  should  have  been  the  same  value  for  the 
modulus  of  elasticity  and  the  Ramberg-Osgood  coefficients  as 
previously  computed.  Since  a  relatively  small  number  of 
points  were  entered  by  hand  (53)  as  compared  to  the  number 
of  points  that  will  be  taken  by  the  program  (1000) ,  larger 
errors  were  expected  on  this  test  than  will  be  expected  on 
an  actual  data  run.  The  results  of  the  test  are  shown  below. 

Kaiser  AXIAL 


E  =  10.117E6 
n  =  21.58 
S  =  1.479343 


E  »  10.113E6 
n  =  23.475 
3  =  1 . 5177E47 


This  large  difference  (especially  in  3)  was  due  to  the  stress- 
strain  pairs  being  calculated  in  stress  steps  of  500  psi. 

The  program,  using  1000  points,  will  have  a  difference  of 
approximately  16  psi  between  points.  This  will  make  the 
accuracy  of  the  program  better  than  the  capability  of  the 
system  to  read  loads  and  strains  accurately. 

The  second  test  of  AXIAL  was  an  actual  specimen  tension 
test.  A  specimen  of  7075-T6  Aluminum  was  fabricated  from 
the  same  sheet  stock  as  the  notched  panel  specimens.  This 
was  inserted  into  the  MTS  machine  which  was  operated  in  a 
constant  strain  rate  mode.  The  specimen  was  loaded  to 
produce  a  strain  of  14,000  microinches/inch,  at  a  strain 
rate  of  1000  microinches/inch/minute.  One  thousand  data 
points  were  taken  at  constant  intervals  of  strain  during 
the  test.  The  program  then  computed  the  modulus  of  elasticity 
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and  the  Ramberg-Osgood  coefficients.  These  were  compared 
to  hand  calculations  made  from  the  Load  vs.  Strain  plot  made 
by  the  MTS.  The  results  of  both  are  listed  below. 

AXIAL  Hand  Calculation 

E  =  9.5042E6  E  =  9.5000E6 

n  =  31.2154  n  =  33.0383 

8  =  3 . 1168E63  3  »  9.7360E66 

A  plot  of  the  actual  data  and  the  points  computed  using  the 
calculated  material  properties  are  shown  in  Figure  6. 

B.  PROGRAM  FLTTST 

The  first  test  performed  with  FLTTST  was  one  to  see  if 
the  random  load  spectrum  followed  and  Mil  Spec  load  spectrum. 
This  was  done  by  making  several  runs  of  100  hour  simulations 
and  comparing  the  spectra  (Table  VI) .  As  would  be  expected, 
the  spectra  did  not  match  exactly  due  to  the  random  way  in 
which  loads  were  chosen.  The  largest  differences  were  in 
the  area  of  the  fewest  loads,  as  also  would  be  expected. 

A  sample  from  a  strip  chart  of  the  Processor  output  and  the 
output  from  a  load  cell  while  the  MTS  is  being  driven  by 
FLTTST  is  shown  in  Figure  7. 

Several  sample  tests  were  run  to  check  the  accuracy  of 
FLTTST  and  DAMAGE.  Knowing  that  a  full  system  test  was 
going  to  be  run,  the  sample  tests  were  run  using  the 
material  properties  of  the  specimen  (Fig.  1)  to  be  used  in 
the  full  scale  tests.  These  material  parameters  and  properties 
are : 
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STRAIN  (;iin/in  x  1000) 
Figure  6.  Axial  Tension  Test 


A  =  . 7136  sq.  in. 

E  =  10.12  E  6  psi 
K  =  2.73 
3*1  479  E  43 
n  »  21.58 
0^  =  72,000  psi 

°ULT  =  82 ' 000  psi 
a_.T  *  60,000  psi 

Jr.u 

TABLE  VI 

SAMPLE  LOAD  SPECTRA 


%  Load 

MIL  SPEC 

100HR 

100HR 

100HR 

100HR 

35 

17000 

1698 

1726 

1662 

1721 

45 

9500 

977 

998 

948 

937 

55 

6500 

658 

637 

690 

646 

65 

4500 

451 

434 

442 

450 

75 

2500 

228 

234 

246 

250 

85 

1360 

133 

139 

146 

131 

95 

440 

49 

21 

43 

53 

105 

150 

13 

8 

15 

8 

115 

40 

5 

1 

6 

3 

125 

16 

0 

2 

2 

1 

42006 

4200 

4200 

4200 

4200 
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igure  7.  Sample  Random  Loading 
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For  these  uses  the  limit  load  was  set  at  the  load  that  would 
cause  local  yield  stress  or  at  two-thirds  of  the  ultimate 
stress,  whichever  was  the  lower.  In  this  case  two-thirds 
the  ultimate  stress  was  the  limiting  factor  which  gave  a 
limit  stress  of  54,667  psi,  produced  by  a  load  of  14,289  lbs. 
Assuming  a  6.5-g  aircraft  limit  load,  this  gave  a  1-g  stress 
of  8409  psi,  produced  by  a  load  of  2198  lbs.  These  are  the 
limits  asked  for  by  the  program  FLTTST  and  all  other  loads 
are  fractions  of  the  limit  load.  The  1-g  load  is  used  as 
the  minimum  load  in  the  fatigue  cycle. 

To  check  the  accuracy  and  speed  of  different  loading 
step  sizes,  sample  calculations  were  performed  and  compared 
with  the  computed  actual  load  that  should  be  applied  (Table  VII) . 
For  this  check  the  material  properties  of  the  sample  specimen 
were  used.  The  conclusions  were  that  a  step  size  of  25 
should  be  used  since  this  gives  good  accuracy  and  takes 
approximately  four  minutes  per  hour  of  flight  time  simulated. 

C.  PROGRAM  DAMAGE 

Further  tests  were  performed  primarily  to  check  the 
DAMAGE  program,  but  the  FLTTST  program  was  used  to  generate 
some  of  the  loads  to  run  through  DAMAGE.  The  first  was  a 
handmade  test  of  ten  cycles  that  had  one  load  that  would 
send  the  specimen  into  the  plastic  region  locally  and  other 
loads  that  would  produce  recordable  damage.  The  purpose  of 
this  check  was  to  check  the  program's  calculation  of  the 
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STEP  SIZE  COMPARISONS 


residual  stresses  and  strains  and  the  difference  in  damage 
due  to  those  residual  stresses  and  strains.  There  were 
three  loads  that  caused  damage  (85%,  95%  and  125%  limit  load) 
and  one  load  that  caused  local  stresses  above  the  propor¬ 
tional  limit  (125%  limit  load) .  Residual  stresses  and 
strains  were  calculated  along  with  the  total  damage  and 
compared  with  the  program's  results  (Table  VIII).  Differences 
in  calculations  were  within  roundoff  error. 


TABLE  VIII 

TEN  CYCLE  LOAD  SPECTRUM 


%  LIMIT  LOAD 

LOAD  (lbs) 

STRESS 

R 

D 

35 

5001.15 

19133 

0.4395 

0 

65 

9387.85 

35532 

0.2367 

0 

45 

6430.05 

24599 

0.3418 

0 

85 

12145.70 

46465 

0.1810 

1.0324E-6 

125 

17861.30 

66944 

0.1048 

3.0193E-5 

35 

5001.15 

17725 

0.3957 

0 

95 

13574.60 

50488 

0.1389 

3 . 8419E-6 

45 

6430.05 

23185 

0.3025 

0 

75 

10716.80 

39567 

0.1772 

0 

55 

7858.95 

28646 

0.2448 

0 

HAND  CALCULATION 

JR  =  1313 

psi 

r> 

II 

7. 3442E-6 

in/in 

D  =  3 . 5067E-5 
PROGRAM  CALCULATION 

=  -1313  psi  e;r  =  7 . 3443E-6  in/in 

D  =  3.5077E-5  7  =  66,944  psi 
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The  last  set  of  tests  ran  the  entire  system.  Program 
FLTTST  was  used  to  drive  the  MTS,  and  output  from  the  MTS 
load  cell  was  read  by  the  program  and  stored  on  tape. 
Photoelastic  pictures  were  taken  of  the  specimen  before  the 
test  and  during  the  test  at  every  load  above  75%  of  the 
limit  load.  A  strip  chart  reading  of  the  program  output 
and  the  load  cell  output  (Fig.  7)  were  read  out  simultan¬ 
eously  to  enable  immediate  analysis  of  the  program  load 
and  MTS  load  correlation. 


TABLE  IX 


THREE  HOUR  FLIGHT  SIMULATION  TEST 


#LOADS 

% LOADS 

STRESS 

R 

D 

2 

85 

46482 

0.1770 

1 . 0726E-6 

3 

95 

52029 

0.1581 

4 . 4607E-6 

2 

105 

57576 

0.1429 

1 . 2222E-5 

HAND 

CALCULATION 

RESULTS 

D  =  3 . 9972E-5 
PROGRAM  RESULTS 
D  =  3 . 9902E-5 


The  first  test  was  a  three  flight  hour  simulation  and 
the  second  was  a  two-hour  simulation.  During  both  tests 
the  highest  load  applied  to  the  specimen  was  105%  of  the 
limit  load.  3y  reference  to  Table  VII,  105%  of  the  limit 


load  does  not  cause  local  stresses  above  the  proportional 
limit;  therefore,  no  residual  stresses  need  to  be  calculated 
and  all  the  1-g  stresses  will  be  the  same.  The  results  of 
the  tests  (Tables  IX  and  X)  show  that  the  program  and  hand 
calculations  agree  within  round-off  error. 

TABLE  X 

TWO  HOUR  FLIGHT  SIMULATION  TEST 


#LOAD 

%LOAD 

STRESS 

R 

D  D 

3 

85 

46482 

0.1770 

1 . 0726E-6 

1 

105 

57576 

0.1429 

1. 2222E-5 

HAND  CALCULATION  RESULTS 
D  =  1. 5440E-5 
PROGRAM  RESULTS 
D  =*  1 . 5534E-5 
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VI.  CONCLUSIONS 


It  is  believed  that  this  system  will  greatly  simplify 
and  increase  the  accuracy  of  predicting  fatigue  damage. 
Almost  all  of  the  calculations  that  would  normally  have  to 
be  done  by  hand,  one  step  at  a  time/  are  now  performed  by  the 
programs.  A  few  of  the  input  values  were  left  as  hand 
calculations  as  a  safety  feature,  so  that  the  user  would  be 
aware  of  the  magnitude  of  the  loads  at  which  the  MTS  would 
be  operating. 

Although  several  tests  and  examples  of  program  verifi¬ 
cation  were  included  in  this  thesis,  many  more  checks  were 
performed  to  insure  that  the  system  was  operating  correctly. 
The  only  area  in  which  problems  occurred  was  with  the 
Program  AXIAL.  If  enough  data  points  are  not  taken,  E  will 
still  be  evaluated  correctly  but  the  Ramberg-Osgood  coeffi¬ 
cients  will  suffer  greatly  in  accuracy. 

Even  though  these  programs  will  greatly  enhance  analysis 
in  the  fatigue  testing  area,  the  potential  of  the  data 
acquisition  system  has  barely  been  tapped.  The  system  as 
it  now  stands  could  handle  over  twice  as  much  data  input 
and  output.  With  added  control  cards,  the  system  could  be 
expanded  at  least  tenfold. 
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APPENDIX  A 
PROGRAM  FLOW  CHARTS 


PROGRAM  AXIAL 
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PROGRAM  FLTTST 


INPUT  CONSTANTS 
AND 

CREATE  DATA  FILE 


SELECT  RANDOM 
NUMBER  SEED 


OUTLOAD=OUTLOAD- 

STEP 


SELECT  NEW 
RANDOM  NUMBER 


SELECT  LOAD 
THAT  CORRESPONDS 
TO  RANDOM  NUMBER 


OUTLOAD*lg  LOAD 


)UTLOAD>SELECTED 
\  LOAD  ^ 


TAKE  DATA  AND 
STORE  ON  TAPE 


TAKE  DATA  AND 
STORE  ON  TAPE 


FINISHED? 


OUTPUT  NUMBER 
OF  EACH  %  LOAD 
SELECTED 
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PROGRAM  DAMAGE 


YES 
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260  PRINT  "  TEST 

270  END 
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PRINT  “THIS  PROGRAM  IS  FOR  AXIALLY  TESTING  A  MATERIAL  SPECIMEN 
PRINT  "  TO  FIND  THE  MODULUS  OF  ELASTICITY,  POISSON"S  RATIO  AND 
PRINT  “THE  RAMSBURG-OSGOOD  COEFFICIENTS  FOR  THE  PLASTIC  REGION 
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PRINT  "THIS  PROGRAM  IS  USED  TO  CALCULATE  FATIGUE  DAMAGE" 
PRINT  "FROM  DATA  OBTAINED  BY  USING  THE  FLIGHT  TEST  PROGRAM. 
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